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Abstract 

Random matrices are used in fields as different as the study of multi-orthogonal 
polynomials or the enumeration of discrete surfaces. Both of them are based on 
the study of a matrix integral. However, this term can be confusing since the 
definition of a matrix integral in these two applications is not the same. These 
two definitions, perturbative and non-perturbative, are discussed in this chapter 
as well as their relation. The so-called loop equations satisfied by integrals over 
random matrices coupled in chain is discussed as well as their recursive solution 
in the perturbative case when the matrices are Hermitean. 



1 Introduction: what is a matrix integral? 

The diversity of aspects of mathematics and physics exposed in the present 
volume witnesses how rich the theory of random matrices can be. This large 
spectrum of applications of random matrices does not only come from the nu- 
merous possible ways to solve it but it is also intrinsically due to the existence, 
and use, of different definitions of the matrix integral giving rise to the partition 
function of the theory under study. 

Back to the original work of Dyson |Dys62| , the study of random matrices 
is aimed at computing integrals over some given set of matrices with respect 
to some probability measure on this set of matrices. In order to be computed, 
these integrals are obviously expected to be convergent. Nevertheless, one of 
the main applications of random matrices in modern physics follows from a 



slightly different definition. Following the work of |Bre78] . the matrix integral 
can be considered, through its expansion around a saddle point of the integrand, 
as a formal power series seen as the generating function of random maps, i.e. 
random surfaces composed of polygons glued by their sides^. Whether this 
formal series has a non- vanishing radius of convergency or not does not make 
any difference: only its coefficients, which take finite values, are meaningful. 

The issue whether these two definitions do coincide or not was not addressed 
for a long time and led to confusions. In particular it led to a puzzling non- 
coincidence of some result in the literature [Akc96bJ IBre99j IKan98j . Their 
computations of the same quantity, even if proved to be right, did not match. 
This puzzled was solved by Bonnet, David and Eynard [BonOO] who were able 
to show that the mismatch between the two solutions is a consequence of the 
discrepancy between the definitions of the matrix integrals taken as partition 
functions. 

Since some of the topics discussed in the present chapter do depend on the 
definition one considers for the partition function whereas some other issues 
do not, section [2] is devoted to the precise definition of these different matrix 
integrals. In section O we present the loop equations which can be used to 
compute the partition function and correlation functions of a large family of 
matrix models. Section [5] is devoted to a review of one the solution of the 
one Hermitean matrix model through the use of the so-called loop equations. 
Section[5]generalizes this method to an arbitrary number of Hermitean matrices 
coupled in chain. Finally, section [6] gives a short overview of generalizations and 
applications of this very universal method. 

2 Convergent vs formal matrix integral 

One of the most interesting features in the study of random matrices is the 
behavior of the statistic of eigenvalues, or correlation functions, as the random 
matrices become arbitrary large. This limit is not only very interesting for 
its applications in physics (study of heavy nuclei, condensed matter...) but 
also in mathematics: the knowledge of the large size limit allows to access the 

1 See chapter 26 for an introduction this topic and |Eyn06] and reference therein for the 
generalization to multi-matrix integrals. 



asymptotics of a large set of multi-orthogonal polynomial^. 

Most of the usual technics used in random matrix theory fail in the study 
of the large matrix limit. However, one possible way to address this problem is 
to try to use naively some saddle point analysis. Let us consider the example of 
a Hermitean one matrix with polynomial potential to illustrate this procedure. 
The partition function is given by the matrix integral: 

Z(V)= [ dMe-T^nM) 

JH.N 

where one integrates over the group Ti^ of Hermitean matrices of size N with 
respect to the measure 

N 

dM : = J] dMu 11 dft (M i:j ) d% (Mij) 

defined as the product of the Lebesgues mesures of the real components of the 

matrix M divided by the volume of the unitary group of size N. For the sake of 

d 

simplicity, one assumes that the potential V{x) = y^^ja; fc+1 is a polynomial. 

k=0 

Notice that the direct saddle point analysis of this integral does not make sense 
in general. 

In order to fix this, let us consider a more general problem. Instead of 
considering Hermitean matrices, we consider normal matrices of size ./V whose 
eigenvalues lie on some arbitrary path 7 in the complex plane: Hn{j) is the 
set of matrices M of size N x N such that there exists U € U(N) and X = 
diag(xi, . . . , ccjv) with Xi € 7 satisfying M = UXU^ . 

With this notation, the set of Hermitean matrices is Hn (M). Given a fixed 
potential V(x), one considers the family of matrix integrals over formal matrices 
on arbitrary contours 7: 

2(V, 7 )= [ dMe-T^vW. 

As in the Hermitean case, one can integrate out the unitary group to turn this 

2 See |Meh04| for a nice review of these application and all the other chapters of the present 
volume. 



partition function into an integral over the eigenvalues of the random matri: 




N 

n -fE y ^) 

Z{V,i)= ... ]Jdx t Y[( Xi - Xj f e i=i 

J f i i=i i<j 

However, given a polynomial potential of degree d + 1, not every path 7 is 
admissible. Indeed, there are only d + 1 directions going to infinity where 
5ft [V(a;)] > as x — > 00 and where the integrand decreases rapidly enough for 
the integral to converge. Thus there exists d homologically independent paths 
on which the integral f dxe~~ v ^ is convergent. Let us choose a basis {^i}f =1 

of such paths. Every admissible path 7 for the eigenvalues of the random matrix 

d 

can thus be decomposed in this basis: 7 = ^^Ki7i- 

i=i 

Using this decomposition, for any admissible path 7, the partition function 
reduces to 

N 

m h} 1 r» = m y: n^T / ni nd no* - x ^ 2 e ,=i 

{ni }i=i nt - ^Ti'x-xt/ i= i iKj 

(2.1) 

where one sums over all integer (i-partitions (ni, . . . , n^) of N, i.e. the sets of d 

integers {rii}f =1 satisfying n\ H + rid = N. 

The requirement of convergence of the integral only fixes the asymptotic 
directions of the paths 7j's. We still have the freedom to choose their behavior 
away from their asymptotic directions. Does there exist one choice better than 
the others? One is interested in performing a saddle point analysis of the matrix 
integral. One thus has to look for the singular points of the action, i.e. the 
solutions of V'{x) = 0. There exist d such solutions i = 1,... ,d, i.e. as 
many as the number of paths 7$ in one basis. In the case of the one matrix 
model with polynomial potential exposed in the present section, it was proved 
following [Ber07] that there exists a good basis in the sense that every path 7, 
is a steepest descent contouiQ. More precisely, along any path 7$, the effective 

3 This procedure can be generalized to multi- matrix models using the HCIZ formula ItzSO 
Har57] presented in chapter 17. 

4 The existence of a good path is conjectured to hold for all other matrix models discussed 
in this chapter. However, the proof is known, at the time these lines are written, only in the 
one matrix model case. 



potential felt by an eigenvalue x, V e ff(x) = V(x) — (In (detx — M)), behaves 
as follows: its real part decreases then stays constant on some interval and then 
increases, whereas its imaginary part is constant then increasing and finally 
constant. 

Such a steepest descent path can thus be seen as a possible vacuum for one 
eigenvalue. Each (/-partition of N hence corresponds to one vacuum for the 
theory, or one saddle configuration for the random matrix. The formula eq l2.ll 
can be understood as a sum over all possible vacua of the theory: 

d K ni 

Z(V,j)= N^-^Z{V,H}\m,...,n d ) 

n±-\ \-n,i=N i=l 

where the partition function with fixed filling fractions = ^ 

N 

, N _ ~fl>(*i) 

Z(V, {7i}K,---,n d ) := / \\dxi T\.( Xi ~ e i=1 

is the weight of a fixed configuration of eigenvalues, or the partition function 
of the theory with a fixed vacuum labeled by a partition (ni, . . . , rid). 

Assuming that the paths ji are good steepest descent paths, the partition 
functions with fixed filling fractions can be computed by saddle point approx- 
imation, i.e. perturbative expansion of the integral around a saddle as t — > 0. 
Further assuming that one can commute the integral and the power series ex- 
pansion, the result is a formal power series in t whose coefficients are gaussian 
matrix integrals: 

Z{V,{li) N>- • • > n rf) ~ 2 formal(V,{li} \ n l: ■ ■ ■ i n d) when t^O 
with 

Z formal := e-^V^±t^fl([ dM ) (jh TrSVm)) 

k=Q 1 K - i=l \ J H ni (7i) J \i=l / 

-i^ v " te)Tr(Ml - ?a '^ )2 
xe i Y[det (Mi ® l n . - l Ui ® Mj) 2 

j<i 

(2.2) 

where {^} i=1 d denote the d solutions of the saddle point equation V'(^i) = 
0, 5Vi(x) := V(x) — V(ii) ^-(x — £j) denotes the non-gaussian part of 



the Taylor expansion of the potential around the saddle and the notation 

d 

II (i/M 7l ) dMi ) stands for the multi P le integral J H ^ M dM x ■ ■ ■ J Hrid ( 7d ) dM d- 

8=1 

This formal series in t is referred to as a formal matrix integral even though it 
is not a matrix integral but a formal power series in t. 

This construction can be thought of as a perturbation theory: the matrix 
integral 2(V) is the non-perturbative partition function of the theory whereas 
the formal matrix integral Zf orma i(V,M\ni, . . . , n<j) is a perturbative partition 
function corresponding to the expansion around a fixed vacuum (m, . . . , rid) m 
the basis (71, ... ,7 d ). 

Since these two possible definitions of the partition function might be con- 
fused, let us emphasize their main differences, concerning their properties as 
well as their applications: 

• The convergent matrix integral is fixed by a choice of potential V together 
with an admissible path 7. The formal matrix integral depends on a po- 
tential V of degree d, a basis of admissible paths {7i}f =1 and a d partition 
of N, Klti- 

• By definition, the non-perturbative partition function is a convergent ma- 
trix integral for arbitrary potential, provided the paths 7, are chosen 
consistently. The perturbative integral is a power series defined for ar- 
bitrary potentials, integration paths and filling fractions. It might be 
non-convergent, and will be for most of combinatorial applications; 

• The logarithm of the perturbative partition function always has a 
expansion, whereas the non-perturative one does not have one most of 
time (see section 13.31) . 



• The formal matrix integral is typically used to solve problems of enumer- 
ative geometry such as enumeration of maps or topological string theory. 
The convergent matrix integral is related for example to the study of 
multi-orthogonal polynomials. 

3 Loop equations 



Even if the perturbative and non-perturbative partition functions do not coin- 
cide in general, they share some common properties. One of the most useful 



is the existence of a set of equations linking the correlation functions of the 



theory: the loop equations. These equations, introduced by Migdal |Mig83 
are simply the Schwinger-Dyson equations applied to the matrix model setup. 
They proved to be an efficient tool for the computation of formal matrix inte- 
grals as the explicit computation of one class of one Hermitean formal matrix 
integral by Ambjorn and al [Amb93| proves. 

3.1 Free energy and correlation functions 

One of the main quantities used in the study of matrix integrals is the free 
energy which is defined as the logarithm of the partition function: 

N 2 

In the formal case, where Z is the generating function of closed discrete surfaces, 
the free energy enumerates only connected such surfaces. 

In order to be able to compute the free energy, but also for their own 
interpretation in combinatorics of maps or string theorjj^l, it is convenient to 
introduce the following correlation functions: 

W k (x 1 ,...,x k ):={Tr 3— Tr 1— . . . Tr ' 



x\ -M x 2 -M x k -M 
where the index c denotes the connected part and 

1 ff (M,-q,f m , 
x-M {x-Ci) k+1 

8=1 k=0 V S V 

has to be understood as a formal power series. It is also useful to introduce the 
polynomial of degree d — 1 in x 

p , x / n^Hw ' i \ 

\ 8=1 I c 

3.2 Loop equations 

The non-perturbative partition function is given by a convergent matrix inte- 
gral. It should thus be invariant under change of the integration variable M 



5 They are generating functions of open surfaces, as opposed to the free energy which 
generates surfaces without boundaries. The interested reader can refer to chapter 31 of the 
present book or the review |Eyn06] for details of this interpretation. 



(or its entries). The name loop equation refers to any equation obtained from 

re partition function under a change 



the invariance to first order in e — ► of t 



of variable of the form M -> M + e8(M 



f dMe-T^V(M) = f d (M + e8(M))e-T'*V(M+ e 6(M)). 

To first order in e, this means that the variation of the action should be com- 
pensated by the Jacobian of the change of variables: 

— < Tr V'(M)5(M)) = (J(M)) . 

Actually, the form of the changes of variable considered is limited to two main 
families of 5{M). This allows to give a recipe to compute the Jacobian rather 
easily as follows. 

• Leibnitz rule: 

J \A{M)B{M)\ = {J [A(M)B(m)}} m _ M + {J [A(m)B(M)}} m ^ M ; 

l-l 

• Split rule: J [A(m) M l B(m)] =J^ Tr [A{ m ) Mj ] Tr [M , - J '- 1 B(m)] ; 

j=0 

i-i 

• Merge rule: J [A(m) Tr (M l B(m))] =J^Tr [A(m)AP B^M 1 -^ 1 ] ; 

i=o 

• if there is no M: J [A(rn)\ = 0. 

The formal matrix integral is obtained from Gaussian convergent integrals 
by algebraic computations which commute with the loop equations thus 

Theorem 3.1 The formal matrix integrals satisfy the same loop equations as 
the convergent matrix integrals. 

3.3 Topological expansion 

The loop equations are a wonderful tool for the study of formal matrix integrals. 
From now on, we restrict our study to these formal power series, leaving aside 
the convergent matrix integrals. 



3 It can be equivalently seen as integration by parts. 



Following an observation originally made by t'Hooft in the study of Feynman 
graphs of QCD [tHo74| . one can see that the exponent of N in the free energy T 
is the Euler characteristic of the surface enumerated by this partition function. 
Thus, T admits a ™ expansion 

oo 
9=0 

commonly called topological expansion since the terms F^ 9 ' of this expansion 
are generating functions of connected closed surfaces of fixed genus g. 

As for the free energy, one can collect together coefficients with the same 
power of N in the correation function and get 

W k (xi,...,x k ) =X)( t) Wj : h) (x l ,...,x k ) 
h=o V ' 



as well as 



P k (x,x x , . . . , x k ) = y~] ( — ) Pp(x,xi, ...,x k ) 

h=0 V 1 J 

where the coefficient are formal power series in t independent of N. 

Both the convergent (non-perturbative) partition function and the formal 
(perturbative) matrix integral are solution to the loop equations. Nevertheless 
they do not coincide in general considered that the loop equations have not a 
unique solution. Indeed, in order to make the solution of these equations unique, 
one has to fix some "initial conditions" satisfied by the sought for solution and 
the convergent and formal matrix integrals are not constrained by the same 
kind of conditions. 

On the one hand, the formal matrix integral has well defined constraints: it 
has a -fij expansion and the small t, large x, limit of any correlation function 
is fixed by the choice of filling fractions. In other words, by fixing the filling 
fractions, one prevents the eigenvalues of the random matrix from tunneling 
from one saddle to another, i.e. from one steepest descent path to another. 
There is no instanton contribution. 

On the other hand, the convergent matrix integral does not admit, in gen- 
eral, any ™ expansion. Moreover, its resolvent is not normalized by any arbi- 
trarily fixed choice of filling fraction: it is rather normalized by some equilibrium 



conditions on the configuration of the eigenvalues, which, thanks to tunneling, 
gives instanton corrections to the classical partition function around the true 
vacuum of the theory. This means that the eigenvalues of the matrix distribute 
on the different paths of the basis in such a way that they are in equilibrium 
under the action of the potential and their mutual logarithmic repulsion. 

In the formal case, one of the main properties of the correlations functions 
is the existence of a topological expansion. Let us plug these topological expan- 
sions into one set of equations obtained by considering the change of variable 
k 

of type SM = x \j Y\ Tr X .^ M ■ They read, order by order in N~ 2 : 
i=i 

V'(x)wi%(x,J) = + 

m=0IcJ 

+ P {h) (x T) + V-3 W n h> ^,J\{^})-Wi h \j) 

in V i ) i / j Q X j x—Xj 
i=l 

(3.2) 

where J stands for {xi, . . . ,x n }. This is the hierarchy of equations which is 
solved in the following section. 

Remark 3.1 Remember that the correlation functions can be seen as the generating 
functions of discrete surfaces of given topology. In this picture, the loop equations get a 
combinatorial interpretation: they summarize all the possible ways of erasing one edge 
from surfaces of a given topology. This gives a recursive relation among generating 
functions of surfaces with different Euler characteristics. This inductive method was 
introduced in the case of triangulated surface by Tutte Tut62J without any matrix 
model's representation of the considered generating functions. 

4 Solution of the loop equations in the 1MM 

The solution of the loop equations in their topological expansion has been 



under intensive study since their introduction by Migdal Mig83 . In particular, 
[Amb93j proposed a general solution of these equations in the one matrix model 
case for the so-called one cut case, i.e. the case where only one of the filling 
fractions e, doesn't vanish. The first steps in the study of the 2-cut case were 
then performed by Akemann in |Ake 96aj . 

Later, in 2004, Eynard |Eyn04| solved the loop equations eq J3.2I for the 
formal integral for an arbitrary number of cuts, i.e. compute all the terms in 
the topological expansion of any correlation function as well as the free energy's 



-^■-expansion for arbitrary filling fractions. This solution relies heavily on the 
existence of an algebraic curve encoding all the properties of the considered 
matrix model: the spectral curve. Let us first remind how the latter can be 
derived. 

4.1 Spectral curve 

Consider eq J3.2I for (h,n) = (0,0): it is a quadratic equation satisfied by the 
genus one point function: 

W ( ®\xf - V'{x)w[°\x) = -P(°)(x) (4.1) 

called the master loop equation. This can be written as an algebraic equation 

H 1MM (xM 0) ) = 

where Himm(x,u) is a polynomial of degree d in x and 2 in y. The algebraic 
equation Himm(%>u) = is the basis of the solution presented in this section 
and will be referred to as the spectral curve of the considered matrix model. 

A first corollary of this equation is the multi-valuedness of W±°\x) as a 
function of x. Indeed, considered P{x) known, one can solve this equation and 
get: 

w (o ){x) _ v'(x)±vv'(xr-4PW(x) _ (42) 

A priori, for any value of the complex variable x, there exist two values of 
Wi°\x). Since, from the definition 14. 2\ its large x behavior is known to be 



d 

£ 



Wl°\x) ~ = ^ as^oo, (4.3) 

one has to select the — sign in order to get the physically meaningful correlation 
function. 

If one can relieve this ambiguity by hand for the genus zero one point func- 
tion, the computation of the complete topological expansion of all the correla- 
tion functions would imply such a choice at each step. 

On the other hand, one can totally get rid of this problem by understanding 
where it originates from. Any correlation function is defined as a formal power 
series both in t — ► and in x — > oo. The coefficients of the -^--expansion of 



Wi(a?) are thus well defined only around x — > oo, as this series might have a 
finite radius of convergency: it is not an analytic function of x. In order to get a 
monovalued function, one has to extend this series further than this radius. The 
master loop equation tells us how one can proceed: instead of considering the 
correlation function W\(x) as a function of the complex variable x, one should 
consider it as a function defined on the spectral curve. That is to say that one 
should not consider W^ l \x) as functions of a complex variable x but rather as 
functions of a complex variable x together with a + or — sign corresponding 
to the choice of one branch of solution of the master loop equation. The tools 
of algebraic geometry are built to be able to deal with such situations and we 
present it in the next section. 

4.2 Algebraic geometry 

Consider an algebraic equation £(x,y) = in y and x of respective degrees 
degree d y + 1 and d x + 1. 

A classical result of algebraic geometry states that there exists a compact 
Riemann surface C and two meromorphic functions x(p) and y(p) on it such 
that: 

Vp G C, £{x(p),y(p)) = 0. 

By abuse of language, one shall use the term spectral curve to denote the 
Riemann surface £, the triple (£, x, y) and the equation £(x, y) = in the 
following, when no ambiguity can occur. 

Let us detail some generalproperties of this spectral curve useful for the 
resolution of the matrix modeQ 

4.2.1 Sheeted structure and branch points 

For a generic fixed value of x, there exists d y + 1 functions of x, y l (x), i = 
0,...,d y solution of the equation £{x,y l {x)) = 0. In other words, a given 
complex number x(p) has d y + 1 preimages p 1 , i = 0, . . . , d y on the surface C 
corresponding to different values of y(p l ). one can thus see the Riemann surface 
C as d y + 1 copies of the Riemann sphere, denoted as x-sheets, glued together, 

7 Most of the properties needed for the study of matrix models can be found in |Far92l 
|Fay73| as well chapter 29 of this volume. 



the function x being injective on each cop; 




How are these sheets glued together to form the Riemann surface £? Two 
sheets merge when two branches of solution in y coincide: y(p l ) = y(p J ) for 
i 7^ j. These critical points aj, called branch points, are characterized by the 
vanishing of the differential dx, i.e. the branch points are solutions of the 
equation dx(ai) = 0. From now on, we suppose that all the branch points are 
simple zeroes of the one form dx. This means that only two sheets merge at 
these points. 

This last assumption implies that, around a branch point a, one has 

y(p) ~ y(a) + y/x(p) - x{a). 

This assumption also implies that, for any branch point a% and any point z 
close to dj, there exists a unique point z such that x(z) = x(z) and z — > as 
z — > ajf). We call z the point conjugated to z around a{. 

The spectral curve C is thus a d y + 1 covering of the Riemann sphere with 
simple ramification points solutions of dx(ai) = 0. 

Example: hyperelliptic curve 

Let us consider an hyperelliptic spectral curve, i.e. a curve given by a 
quadratic equation d y = 1, as in the one Hermitean matrix model case: 

d 

£(x,y) = V 2 - IJO - x(ai))(x - x(bi)) = y 2 - a(x) 
i=i 

where d x := 2d to match the notations of the previous section. The correspond- 
ing Riemann surface can be seen as a two sheeted cover of the Riemann sphere: 
one sheet corresponding to the branch y(x) = Wa{x), and the other one to the 
other branch y(x) = —^Ja{x). These two sheets merge when y(x) takes the 
same value on both sheets, i.e. when y{x) vanishes. The branch points are thus 
the preimages of the points x(at) and x(bi) on the spectral curve. The latter 
can thus be described as two copies of CP 1 glued by d cuts [a^, fej]. 

8 Each copy of the Riemann sphere corresponds to a branch of solution in y of the equation 
S(x,y) = 

9 The application z — » ~z is defined only locally around the branch points and depends on 
the branch point considered and the notation z is abusive. Nevertheless, this application will 
always be used in the vicinity of a branch point in such a way that no ambiguity will occur. 



4.2.2 Genus and cycles 

Generically, the compact Riemann surface C associated to an algebraic equa- 
tion may have non vanishing genus g, and it will be the case in most of the 
applications of the present chapter. 

The Riemann-Hurwitz theorem allows us to get this genus out of the branched 
covering picture of the spectral curve. For example, if there are only simple 
ramification points, it states that 

number of branch points 
9 = -d y + - 

where d y + 1 is the number of x-sheets. 

If the Riemann surface has non-vanishing genus, i.e. it is not conformally 
equivalent to CP 1 , there exist non-contractible cycles on it. In order to deal 
with them, it will be useful to choose a canonical homology basis of cycles 
(Al, . . . , A g , B±, . . . , Bg) satisfying the intersection conditions 

Vi,j = l...,g, Aif)Ai = Btf)Bj = , Aif]B j = S iJ . 
Example: hyperelliptic curve 

Let us keep on considering the example of an hperelliptic curve with 2d 
branch points. From the Riemann-Hurwitz theorem, it has genus g = d — 1. 
This follows the intuitive picture of two Riemann sphere glued by d segments 
giving rise to a genus d — 1 surface. 

One can also explicit a canonical homology basis as follows. First choose 
one cut, for example [ai,&i] and one sheet, e.g. the sheet corresponding to 
the minus sign of the square root. Then define the A-cycle as the cycles on 
the chosen sheet around the cut [aj+i , • Finally, define the £?j-cycle as the 
composition of the segment [a\, aj-fi] in the chosen sheet and [oj+i, a±] followed 
in the opposite direction in the other sheet (see fig fT] for the simplest example 
of a genus 1 surface). 

4.2.3 Differentials 

The meromorphic differentials on the Riemann surface L and their properties 
will play a crucial role in the following. In particular, let us remind a fundamen- 
tal result concerning meromorphic differentials: a meromorphic differential df 
on Riemann surface C of genus g equipped with a basis of cycles {Ai, Z3j}f =1 ,is 




Figure 1: Genus 1 hyperelliptic curve: it is built as two copies of the Riemann 
sphere glued by two cuts [01,02] and [03,04]. The unique ^4-cycle encircles 
[03,04] while the £>-cycles goes through both cuts. 

defined uniquely by its ,4-cycles f A df on the one hand and its singular behav- 
ior, i.e. the position of its poles and the divergent part of its Laurent expansion 
around the latter. 

For example, one introduce one of the main character of the resolution of 
loop equations as follows: 

Definition 4.1 Let the Bergman kernel B(p,q) be the unique bi- differential in 
p and q on C defined by the constraints as a differential in p: 

• it has a unique pole located at p — > q which is double without residue. In 
local coordinates, it reads 

B(p, q) ~ -, — + regular when p — > q; 
(p -qY 

• it has vanishing A-cycle integrals: 

V* = l,. ..,<?, / B(p,q)=0. 

It is also useful to define the primitive of the Bergman kernel: 

r-P2 

dS PuP2 (q) = / B(z,q) 
Jz=pi 

which is a one form in q with simple poles in q — > p\ and q — > P2 with respective 
residues —1 and +1. 



4.3 The one point function and the spectral curve 

With these few elements of algebraic geometry in hand, let us complete our 
study of the spectral curve of the Hermitean one matrix model. Up to now, one 
has obtained that W^\x) is solution of a quadratic equation which depends 
on a polynomial P(°\x) of degree d — 1, i.e. d variables remain to be fixed. 



From the definition eq l3.ll of the corre 
cycles as circle, independent of t around £. 



ation functions, considering the A%- 
one gets d constraints 



Vi = l,. ..,(*, W [ Q \x)dx= T ^r 



2ivr J Al N 

allowing to fix the coefficients of the polynomial P^(x), since the contour 
integral £4 pics only one residue at £j. One thus gets all the parameters of the 
spectral curve as well as the one point function W± (x). 

Properties of the one matrix model's spectral curve 

The polynomial H\mm{x, y) has degree 2 in y. This means that the embed- 
ding of Cimm is composed by 2 copies of the Riemann sphere glued by g + 1 
cuts so that the resulting Riemann surface Cimm has genus g. Each copy of 
the Riemann sphere corresponds to one particular branch of the solutions of 
the equation Hxmm(,%,v) = 0. Since there are only two sheets in involution, 
this spectral curve is said to be hyper elliptic. It also means that the application 
z — ► z is globally defined since it is the map which exchanges both sheets, i.e. 
which exchange the two branches of the square root in 14.21 

The Riemann surface Cimm has genus g lower than d — 10. 

The function x{z) on the Riemann surface C\mm has two simple poles (call 
them a + and a_), one in each sheet. Near a±, y(z) behaves like: 

y{z) ~ J- +0 (l/x(zf) 

z^a + x(Z) 

and 

y(z) ~ t d x d {z) + 0(x d -\z)). 



10 Indeed, when t — ► 0, the cuts are reduced to double points at £;'s. As t grows, these 
double points give rise to cuts of length of order 2i_ . 

11 Notice that d — 1 is an upper bound. It might happen that two branch points coincide 
resulting in the closing of one cut and decreasing of the genus by one. For some special value 
of the coefficients of the polynomial H\mm one can even get a genus zero spectral curve. For 
application of matrix models to enumeration of surfaces, this very non-generic constraint is 
almost always satisfied (see 'Ey n06| for further considerations on this point). 



4.4 Two point function 

Let us go one step further and consider the loop equation (|3.2p for k = 2 and 
h = 0. It allows to obtain a formula for W^ix, x\): 

a wW(x)-wW(»i) p(Q)/,. T n 
2{V'{x)-W[°\x)) 

A first look at this expression allows to see that this function is multivalued in 
term of the complex variable x and x\. However, one can lift it to a monovalued 
function, actually a 2-form, on the spectral curve by defining 

uj^\ z -, z i) '■= W^\x{z), x{z\))dx(z)dx(zi). 

Q^oj(z, z\) is thus a meromorphic bi-differential on C. One can then study all 
possible singularities of this formula and see that Q^\z, z\) has poles only at 
z — ► z\. On the other hand, the normalization of the two point function around 
the *4-cycles reads §, uj^\z, z\) = for i = 1, . . . , d. These two conditions 
imply that Q^\z,zi) is given by the Bergman kernel (see for example section 
5.2.3 of |Kyn09| ) 

~(o), \ D , _ x „, x dx(z)dx{zi) 
4 (z,z 1 ) = -B(z,z l ) = B(z,z 1 )- {x{z) _ x{zi)f . 

4.5 Correlation functions 

We have now everything in hand to compute any correlation function by solving 
the loop equations. First of all, the study of the one and two point functions 
proved that it is more convenient to promote the multivalued functions Wn 
on the complex plane to monovalued meromorphic forms on 



12. 



^ ._uK*)^ ^fl^M-LX.. n x. „ dx( Zl )dx(z 2 ) 



wHy[zi, . . . ,z n ) := W^ h \zi, . . . ,z n ) Y\_dx(zi) + 5 n ,z&h,Q 

and 



y(z)dx(z) := w[ Q \z)dz 



. It is important to remember that the physical quantities encoded in the 
correlation functions are obtained as the terms of the expansion of the latter 
when their variables approach the physical pole a+ of the spectral curve. 



12 The monovaluedness of the differential form uj& on the spectral curve is obtained by 



induction on the Euler characteristic 2h + n — 2 through the use of the loop equations (|3.2 



Prom the loop equations (|3.2p . one can prove by induction that Un(z\, . . . , z n ) 
with 2h + n > 3 can have pole neither at coinciding points x{zi) = x(zj), neither 
at the poles of x nor at the double points. It may have poles only at the branch 
points. 

Let us now write down the Cauchy formula on the spectral curve: 
wf+i( z . z i'---)4) = Res dS z ' j0 (z)uj^l 1 (z',Zi,...,z n ) 

z'^z 

where o is an arbitrary point of C. Since ui^h^z', z\, . . . , z n ) has poles only 
at the branch point aj, moving the integration contours on C (and not C!), 
one gets contributions from the latter and the boundaries of the fundamental 
domain of £ according to Riemann bilinear formula [Far92] 



vlL(z,Zi,. ..,z n ) = -V" Res dSz^otz)^)^ 1 ,z\, ... ,z n ) 

'—' Z'—Hli 
I 

9 

+E [i'eAi B ( z i z ') i'eBi 4+1 «i» ■ ■ ■ ,*n) 

i=l 

+ b'eBi B ( Z ' Z ') §z'&Ai u n+l ( Z 'i Z ^---i Z n) 

Since the correlation functions and the Bergmann kernel have vanishing *4-cycle 
integrals, the second and third line vanish. One can then plug the expression 
for u^^z', zi, ■ ■ ■ , z n ) coming from the loop equations f|3.2|) into this formula. 
Since the polynomial pjf^ 1 (x(z'), Z\,..., z n ) is regular at the branch points, it 
does not give any contribution and one gets the recursion formula 

^n+l(z,zi,...,z n ) = ^Res z >^ ai K{z,z')[uJ { ^~2 ) {z l ,z',zi,...,z n ) 

i 

h 

+E E u m+iV> ^hShG 7 . • • • ' Zn} 

j=0 lc{zi,...,z„} 

where the sign ^ means that the sum does not involve the terms with (j, |/|) = 
(0,0) or (j, |/|) = (h,n) and the recursion kernel is 

dS , -f(z) 

K(z,z')-, 



2{y(z>) - y{z'))dx{z') 

It is easy to see that this recursive relation on 2h + n — 2, i.e. provided that 
lo^\z, z\) = B(z, z{) is known, it determines all the other correlation functions 
through their topological expansion. 



Remark 4.1 This recursion can be graphically represented in such a way that it 
becomes very easy to remember and allows to recover some of the properties of the 
correlation functions using only diagrammatic proofs. Details on this diagrammatic 
representation can be found in |Eyn09| . 

4.6 Free energies 

In the preceding section, we have been able to compute the topological expan- 
sion of any correlation function W n for n > 0. Let us now address the case 
n = 0, that is to say the computation of the topological expansion of the free 
energy. 

For this purpose, one can build an operator acting from the space of n + 1- 
differentials on C into the space of n-differentials mapping the n + 1-point 
function to the n-point function 

Theorem 4.1 For any h and n satisfying 2 — 2h — n < and any primitive <3> 
of ydx, one has 

w$i\zi, ...,z n ) = 2 _2 h _ n £fl *(*) w 2i( z > *i» " " " ' Zn )' 

i 

One can guess that this definition can be extended to n = in order to get 
the topological expansion of the free energies as follows: 

Theorem 4.2 The terms of the topological expansion of the free energy of the 
Hermitean one matrix model are given by: 

i 

for h>2. 

This guess can be proved to be right by looking at the derivative of the result 
with respect to all the moduli of the formal integral, i.e. the coefficient of the 
potential and the filling fractions. Indeed, they match the expected variations 
of the free energies when varying these moduli [Che0 6aj. 

5 Matrices coupled in a chain plus external field 

It is remarkable that the recursive formula giving the topological expansion 
of the free energy and the correlation functions depends on the moduli of the 
model only through the spectral curve. One can thus wonder wether the same 



procedure can be applied to solve other matrix models which are known to be 
related to a spectral curve. This is indeed the case for the model of two matrices 
coupled in chain }Che06b] but also for the an arbitrary long chain of matrices 
in an external field |Eyn08|. 

In order to deal with a large family of Hermitean matrix models at once, 
let us consider an arbitrarily long sequence of matrices coupled in chain and 
submitted to the action of an external field. 

The partition function is given by the chain of matrices formal matrix inte- 
gral: 



■Tr y]v k (M fc )~y^c fc , fc+1 M k M k+l 
^chain = / r Vfc=l *=1 / dM, ... dM„ 



t 



formal 

where the integral is a formal integral in the sense of the preceding sectional 
M m+ i is a constant given diagonal matrix M m+ i = A with s distinct eigenvalues 
Aj with multiplicities If. 

(h h Is 

Ai, . . . , Ai, . . . , Aj, . . . , Aj, . . . , A s , . . . , A s 

with = N and one considers the m polynomial potential:] 

i 

d k +l 

V k (x) = -J2 tj fx j - 
i=2 

As in the one matrix model case, the definition of the formal integral requires 
to choose around which saddle point one expands. Saddle points are solutions 
of the set of equations 

V& = 1, . . . , m, Vfc(£ fc ) = c k ^ 1;k Ck-i + c k) k+i£k+i , 3j, ^ m+1 = Xj 

which can be reduced to an algebraic equation with D = sdidz ■ ■ ■ d m solutions. 

This choice is thus equivalent to the choice of a D-partition (m, . . . , no) of 
N giving rise to the filling fractions: 

i%i 

€j = t 

N 



13 The formal integral is a power series in t whose coefficients are Gaussian integrals. See 
Eyn06 for a review on this topic 

It is possible to generalize all this section to potentials whose derivative are arbitrary 
rational functions without any significant modification of the present procedure. 



for %= 1,...,D with D = d\d% . . . d m s and nt arbitrary integers satisfying 

i 

5.0.1 Definition of the correlation functions 

The loop equations of the chain of matrices were derived in |Eyn03 Eyn08| , 
and they require the definition of several quantities. 

oo 

For convenience, we introduce Gj(xj) := X ._ M . = / .~kfr as a formal power 

k=o i 

series in X{ — » oo as well as a polynomial in x, Q(x) = — - — ^ a > where 

S 

S(x) is the minimal polynomial of A, S(x) = Y\_( x ~ ^ e a ^ so define the 

i=l 

polynomials fa(xi, ... , xj) by fa = if j < i - 1, fa-i = 1, and 

/ V({xi) -c iti+ ix i+1 \ 

Ji,j(p^i: ■ ■ ■ > %j ) — Get 

— Cj-ljXj 

\ -Cj-ijXj-i V!(xj) j 

if j > i. The latter satisfy the recursion relations 



i—\,ifi,j ip^i j • • • > -Ej) — i^i) fi+l,j > • • • ; ) 2>t+l fi+2 j • • • j )• 



Let us finally define the correlation functions and some useful auxiliary 
functions. In the following PoLj/(x) refers to the polynomial part of f(x) as 
x — > oo. For i = 2, . . . , m, we define 

Wi^l , Xj , . . . , X m , Z) . — 

) ( Tr (Gi(si)Gi(xi) . . . G m {x m )Q{z))) , 
which is a polynomial in variables Xi, . . . , x m , z, but not in x\, for i = 1, 
VFi(a;i,X2,...,2; m ,2:) := 

Polxi,...,a: m /i,m(a;i, • • • , a? m ) ( Tr {G 1 (x 1 )G 2 (x 2 ) ■ ■ ■ G m (x m )Q(z))) 
which is a polynomial in all variables and, for i = 0, Wo(x) = (Tr G\(x)). We 
also define: 

Xi, . . . ,X mi Z\ X-y) '.= 

PoL ej) ..., ! r m /i,m(»i> • • • ,x m ) (Tr (Gi^)) Tr (Gi(xi)Gj(xj) . . . G m {x m )Q(z))) c . 
All these functions admit a topological expansion: 

Wi = YsiN/t) 1 ' 29 ^ and = ^ N ^~ 2gw S- 



5.0.2 Loop equations and spectral curve 

In this model, the master loop equation reads |Eyn03 Eyn08| : 



W 2 -,i(xi, . . .,x m+ i;xi) + jfWiixi, . . . ,x m+ i) - (V{(xi) - c 1;2 x 2 ) S(x m+ i) 
+(c 1>2 x 2 - V{(x\) + j?W (xi))(j r W 2 (xi,...,x m+1 ) - S(x m+1 )\ 

in 

= i 7/. ( v i( x i) - <H-l,iXi-i ~ c iji+ ix i+ i) Wi+i(xi,Xi, . . . ,x m+ i). 

i=2 

(5.1) 

Let us consider specific values for the variables Xi in order to turn it into an 
equation involving only x\ and x 2 . One defines {xi(x\, x^}™^ 1 as functions of 
the two first variables x\ and x 2 by 

c iti+1 x l+ i(xi,x 2 ) = V-(xi(xi,x 2 )) - a-i,iXi-i(xi,x 2 ). (5.2) 

for i = 2, . . . ,m with the initial conditions x%{x\, x 2 ) = x\ and x 2 {x\, x 2 ) = x 2 . 

Choosing Xi = X{(x\,x 2 ), reduces the master loop equation to an equation 
in x\ and x 2 : 

W / 2 ;i(xi,x 2 ;xi) + — {c x , 2 x 2 - Y(xi)) U(xi,x 2 ) = E(xi,x 2 ) 

where Y(x) = V{{x) - jjW (x), the hat means that the functions are considered 
at the value X{ = Xi(xi,x 2 ), i.e. (f)(xi,x 2 ) := f(xi,x 2 ,X3,X4, . . . , x n ) for an 
arbitrary function /, and one has defined 

N 

U(xi,x 2 ) = W 2 (x 1 ,x 2 ,x 3 , . . .,x m+1 ) - —S(x m+ i), 

and 

E(xi,x 2 ) = -—Wi(xi,x 2 ) + (V{(xi) - c lj2 x 2 ) S(xi,x 2 ). 
Finally, the leading order in the topological expansion gives 

,x 2 ) = (c\ y2 x 2 -Y^(xi)^j U^(xi,x 2 ) (5.3) 

where one should notice that Wi(xi,x 2 ), and thus E(xi,x 2 ), is a polynomial 
in both x\ and x 2 . 

Again, this equation is valid for any x\ and x 2 , and, if we choose x 2 such 
that ci )2 x 2 = Y^(xi), we get: 

H chain ( Xl ,x 2 ) := E<®(x 1 ,x 2 ) = 0. (5.4) 



This algebraic equation is the spectral curve of our model. 



Study of the spectral curve 

The algebraic plane curve H c h a i n {x\,x 2 ) = 0, can be parameterized by a 
variable z living on a compact Riemann surface C c h a in of some genus g, and 
two meromorphic functions x\{z) and x 2 (z) on it. Let us study it in greater 
details. 

The polynomial H chain (xi ,x 2 ) has degree (resp. di + ^j) in x 2 (resp. 

x{). This means that the embedding of £ c hain is composed by 1 + ^ (resp. 

copies of the Riemann sphere, called xi-sheets (resp. X2-sheets), glued 
by cuts so that the resulting Riemann surface C c h a in has genus g. Each copy 
of the Riemann sphere corresponds to one particular branch of the solutions of 
the equation H c h a in(xi, x 2 ) = in x 2 (resp. x\). 

The Riemann surface Cchain has genus g lower than D — s with D = 
sd\ . ..d m . 

One can consider all the variables Xi(p) "■= Xi(x\{p), x 2 (p)) as meromorphic 
functions on £ c h a in as opposed to only x and y in the one matrix model case. 
Their negative divisors are given by 

s 

[xkip)}- = -rkoo - s k ^ Xj 

i=l 

where oo is the only point of C c h a in where x\ has a simple pole, the \\ are the 
preimages of Aj under the map x m +i(p), x m +i(Xi) = Aj, and the degrees r k and 
s k are integers given by r% := 1, r k := d\d 2 . . . d k -i, s m+ i := 0, s m := 1, and 
Sk '■= d k+ \d k+2 . . .d m s. 

Note that the presence of an external matrix creates as many poles as the 
number of distinct eigenvalues of this external matrix M m+ i = 

5.0.3 Solution of the loop equations 

The procedure used to solve the loop equations in the one matrix model cannot 
be generalized in this setup, mainly because the involution z — > z is not globally 
defined on the spectral curve. However, the loop equations can been solved 
using a detour |Eyn08| . This resolution proceeds in three steps. One first shows 

15 The cases of matrix models without external field correspond to a totally degenerate 
external matrix A = eld with only 1-eigenvalue. There are thus two poles as in the 1 matrix 
model studied earlier. 



that the loop equations eq l5.ll have a unique solution admitting a topological 
expansion. One then finds a solution of these equations: 

Lemma 5.1 

E(x(z),y) = —K n (j[(y - V[{x{z*)) + 1 Tr " (5.5) 

where K is a constant and the inverted comas " < . > " means that every time 



one encounters a two point function ^ Tr x ^_ Ml 1^ X ( z j)-Mi / > one re P^ aces 
it by W$(zi; Zj) := ( Tr x{z ^_ Mi Tr x{z3 ]_ Mi + {x{Zi) l x ( Z]) y± ■ 

The matching of the coefficients of the polynomials in y in the left- and right 
hand sides of 15.51 and a few algebro-geometrical computations allows to solve 
the loop equationj^l to get 

Theorem 5.1 The correlation functions of a chain of matrices formal integral 
are recursively obtained by computing residues on £ c hain : 

J^l 1 {z,z 1 ,...,z n ) = ^Res z ,^ ai K(z,z')luj { ^~2 ) {z\z',zi,...,z n ) 

i 

h 

j=0 IC{zi, ...,z„} 

(5.6) 

where, as in the preceding section, 



wWfzi,...,^) = Res N n+2h - 3 { IT TrGi(xi) ) tto(zi) . . . dx(z n 

\i=l / c 

c?x(2:i)(ix(z2) 



(x(zi) - x(z 2 )) 2 

i/ie too poini function uj^ is the Bergman kernel of the spectral curve C c hain, 
the recursion kernel is 

K(z,z ) := -= , 

2(y(z>) - y(z>))dx(z>y 

x and y are two meromorphic functions on C c h a in such that H c f ia i n (x(z),y(z)) = 
for any point z £ C c hain an d a>% are the x-branch points, i.e. solutions to 
dx{ai) = 0. 



5 See [Che06bl |Eyn08l for the detailed proof. 



Remember that H c h a in was defined in 15.41 by H c h a in{xi-, C12X2) = for C12X2 = 
Y(°'(xi). Thus the function x and y can be thought of as continuation to the 
whole spectral curve of x\ and C12X2 respectively. 

The free energy is also obtained by using the same formula as in the one 
matrix case: 

Theorem 5.2 For h > 1 and any primitive <5 of ydx, one has 

i 

Thus, the solution of any chain- matrix model with an external field is ob- 
tained by the exact same formula as the solution of the one matrix model: the 
only difference is the spectral curve used to apply this recursion. 

6 Generalization: topological recursion 

We have seen that the loop equation method gives a unique solution for a large 
family of formal matrix models. The only input of this solution is the spectral 
curve of the considered model. In |Eyn07|, it has been proposed to use equations 
15.61 and l5.7l to associate infinite sets of correlation functions and free energies to 
any spectral curve (£, x, y) where £ is a compact Riemann surface and x and 
y two functions analytic in some open domain of C. 

The free energies and correlation functions built from this recursive pro- 
cedure show many interesting properties such as invariance under a large set 
of transformations of the spectral curve, special geometry relations, modular 
invariance or integrable properties. In particular, it is a very convenient tool to 
study critical regimes and get the universal properties of the matrix integrals 
described in the chapter 6 of this volume. It is also very useful to compare 
different matrix integrals. Eventually, this procedure proved to be efficient in 
the resolution of many problems of enumerative geometry or statistical physics 
such as string theory, Gromov-Witten invariants theory, Hurwitz theory or ex- 
clusion processes such as TASEP or PASEP. Most of the results proposed by 
this approach are still conjectures up to now but the numerous checks passed 
so far tends to prove that this generalization of the loop equation method is a 
very promising 

fielfl 

17 For a review on this subject, see |Eyn09| and references therein. 



The inductive procedure presented in this chapter only allows to compute 
one particular set of observables of multi-matrix models. It does not compute 
correlation functions involving more than one type of matrix inside the same 
trace. These more complicated objects are very important for their application 
to quantum gravity or conformal field theories where they correspond to the 
insertion of boundary operators. In the two matrix model, the loop equation 
method allowed to compute the topological expansion of any of these operators 
|Eyn08| . In the chain of matrices case, only a few of them were computed in 
their large N limit only |Ey n03|. The computation of any observable of the 
chain of matrices is still an open problem which is very likely to be solved by 
the use of loop equations. 

Acknowledgements: It is a pleasure to thank Bertrand Eynard who devel- 
oped most of the material exposed in this chapter and patiently taught me all 
I know about these topics. 
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